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EELATIONS BETWEEN THE ECCENTRICITIES AND INCLI¬ 
NATIONS OF THE ORBITS OF THE PLANETS JUPITER 
AND SATURN. 

By R. T. A. Innes, F.R.S.S.Af. 

(Read June 21, 1911.) 

It is well known that the eight major planets revolve around the Sun 
in nearly the same plane and in elliptical orbits which do not depart 
greatly from circularity. 

In the cases of the Earth, Mars, Jupiter, and Saturn, the equators of 
these planets are not greatly inclined to the common plane of the orbits; 
the equators of Mercury, Yenus, Uranus, and Neptune have not been 
seen, but from the motion of the satellites of Uranus and Neptune it is 
inferred that their equators have large inclinations. 

The precise calculations of astronomers, which unfortunately can only 
extend to a few hundreds of years, agree with observation in showing that 
the eccentricities and inclinations of the various orbits are changing 
slowly. Lagrange, Laplace, and Poisson have proved that the mean- 
distances of the planets are essentially invariable, and the two former 
have, by very rough and inexact methods, shown that the sums of simple 
functions of the eccentricities and inclinations of all the planets will 
always remain small, but on account of the preponderating masses of the 
four outer planets the equations prove nothing as to the stability of the 
orbits of the four inner planets, of which group the Earth is one. It> 
cannot to-day be proved that the orbit of such a planet as the Earth is 
stable. The large uniformity of flora and fauna for great ages gives rise to 
a no doubt well-founded belief that the Earth’s orbit is stable, although 
there may be no mathematical proof yet available. 

The case is different with the planets Jupiter and Saturn ; the pertur¬ 
bations which these planets undergo, due to the presence of the other 
planets, are so insignificant compared to their own action on each other 
that they can be considered apart. In fact, the preponderating masses of 
the Sun, Jupiter, and Saturn in the solar system make these three bodies 
a nearly ideal case of the “ problem of three bodies,” and in consequence 
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of this the rigorous integrals given by Jacobi in his celebrated paper on 
the “ Elimination of the Nodes ” are closely approximated to.* 

If we adopt the notation for Jupiter— 

t = inclination of orbit to the invariable plane, 
e — eccentricity of orbit and S <p = e, 
ra = mass, 

a = ratio of Jupiter’s mean distance to Saturn’s, 


and use accented symbols for the same quantities when applied to Saturn, 
and take the constant of attraction as unity, we have (see Charlier, 
“Mechanik des Himmels ” formulae of 14,* p. 274 of vol. i.)— 


in which— 


/3 \!aG(j)Si — l3 I G(p I St 1 = 0 
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and C and S are written for cosine and sine, and c 1 is a constant. 
Taking— 

1/m = 1047-3, 1/7^ = 3501*6, 

Lx/a = 9-8682770-10, 

and dividing the equation (A) by j3 J a, we get— 

CtySt - qCfaStt = 0 ~Lq = 9-6076936-10 
G<pGc + qG^Gi, = c q = 0-4052225, 


in which both q and c are absolute constants. 

Dr. G. W. Hill found by a comparison between the places given by 
his theory of Jupiter and Saturn and those actually observed that the 
reciprocals of the masses were 1047-38 and 3502-2 respectively; these 
values would change the value of L q to 9-6076314-10, so that one must 
not lay too much stress on the value of this constant. 

If we square and add the two equations, we get— 


G 2 (f> -j- 2gC0C0 x C(i -j- f-i) q 2 G 2 (j> I — c 2 . 


* The not too great inequality of the three masses is also essential, and exists in the 
case of these three bodies—the integrals would furnish no useful information in the case 
of Jupiter, Sun, and Jupiter’s 8th satellite, because of the smallness of the mass of JVIII, 
compared with Jupiter and the Sun. 




Relations between Jupiter and Saturn. 


289 


We take the following figures from Hill’s tables :— 

Le = 8-6835317-10 LC0 = 9-9994938-10 

Le z = 8-7486589-10 LCfr = 9-9993165-10 

t + tl = l° 15' 20-9" LC(t + i z ) =9-9998957-10 

and thus have for the value of the constant c — 

Lc = 0-1471664 c = 1*4033512 
c 2 = 1-9693950 


We can now make a rough discussion of the last equation, viz.— 

00 + 0-81 C0C0 I C(i + i 1 ) + O-16 O0, = 1-969. 

As no cosine can exceed unity, and C0, C(f> zt and C(« + i z ) are essentially 
positive, it is at once evident that at all times the C quantities must each 
be less than but yet exceedingly close to unity. If each is taken as unity 
the equation fails to balance, but only just fails, for we then have— 

1 + 0-810445 + 0-164205 -1*974650, 
whilst c 2 is actually =1-969395. 


Thus the eccentricities and inclinations of the orbits of these two 
planets can never deviate greatly from their present small values. 

Both i and t z can be determined in terms of e and e z , and vice versa. 
We have, in fact— 


^ faCVi + C0 —c) (qG(f) 1 — C0 + c ) 


Six/2 


-V m 


(qO<p z + G<p-c) (0<j> — qG<t> z + o) 


icqG(f> l 


°^ _C S (, + „) 

n c St 

C ^~q S(T+ tl ) 


The first pair of these equations are ill-suited for finding the inclina¬ 
tions, because gC0 r + C0 is very nearly equal to c. It is better to find 
i and t z by trial from the first equation of (A) ; we thus find— 


i = 21' 16-6" 

t z = 54' 4-3" 
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If we adopt the variations in e and e x and i and i x found by other 
investigators,* we can draw some further inference from the equation (A). 
Hill (“Collected Works,” iv., p. 127) has discussed the maximum and 
minimum values of e and e xi and he points out that the maximum eccen¬ 
tricity of one orbit corresponds with the minimum eccentricity of the 
other, and vice versa. From Hill’s figures we get— 

LC<£ (min.)-9-9991924-10 LC0 x (max.) = 9-9999705-10 
LC</> (max.) = 9*9998597-10 LC^ (min.) = 9-9984324-10 

If we write the first equation of (A) thus— 

and substitute the above values of C <j> and C<p iy we get— 

Si = 0-4044 Si r . Min. of variable factor. 

= 0-4051 S ij. Present value (1850). 

= 0-4054 Si x . Max. value. 

These equations show that i and i x increase or decrease together, and 
that with considerable precision we have at all times— 

Ai = 0*405Ai i . 

It must be remembered that Ai and Ai x are the variations of the 
inclinations to the invariable plane of the system. 

Stockwell (“ Smithsonian Contributions to Knowledge,” xviii., 1870) 
gives the following figures :— 

Maximum value of i + i r = 1° 29' 35" 

Minimum value of i + h = 1 1 39 

to which we add the present value =1 15 20*9 (1850) 

Value 2,000 years hence (Le Verrier) ==1 15 1 (3850) 

These values show that with great accuracy we can write the second 
equation of (A) in the approximate form— 

„ _ j (7 = 0-405 

Cf+qCq^ c } c = 1-403 

or making e and e x variable— 

eAe + ^jAe^O. 

* The argument is complete without these, and logically they should not be used ; but 
their use simplifies the reasoning. 


Relations betiueen Jupiter and Saturn. 
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This simple formula shows, as Hill has already pointed out from other 
considerations, that as e increases or diminishes e x decreases or increases, 
and that when e is at a maximum e 1 is at a minimum, and vice versa. It 
further shows that the variations of the eccentricities preserve the ratio 

— 0*405 x^. Hill’s Tables of Jupiter, and Saturn do not give the varia¬ 
tions of the eccentricities directly, but from them I have found the 
centennial variations to be— 

Jupiter +34*06" 

Saturn - 71*63 

These figures give a ratio of — 0*409^, which is a close accordance,, 

especially when it is remembered that Hill’s variations include the action 
of Uranus and the other planets. 

It is remarkable that the ratio of the mean motion of Jupiter to that 
of Saturn ( = 0*4028) is so close to the value of q ( = 0*4052) and that the 
present ratio of the great inequality in the mean longitude of Saturn fo> 
the corresponding inequality of Jupiter (by Hill’s tables =0*4113) is again 
not widely different from q. These approximations may be fortuitous. 


Summary. 

The following relations, which are practically rigorous, connect the 
variations of the eccentricities and inclinations to the invariable plane of 
the planets Jupiter and Saturn :— 

Jupiter. Saturn. 

Ae = -0-405-Ae t . 

e 

A t = +0*405A tl . 
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